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Laser-induced cooling of broadband heat reservoirs 
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We explore, theoretically and experimentally, a method for cooling a broadband heat reservoir, via 
its laser-assisted collisions with two-level atoms followed by their fluorescence. This method is shown 
to be advantageous compared to existing laser-cooling methods in terms of its cooling efficiency, the 
lowest attainable temperature for broadband baths and its versatility: it can cool down any heat 
reservoir, provided the laser is red-detuned from the atomic resonance. It is applicable to cooling 
down both dense gaseous and condensed media. 
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Introduction The advancement of quantum technolo¬ 
gies is demanding new coolingmethods as a follow¬ 
up on existing laser cooling [sHol. A cooling method is 
assessed by its cooling power, thermodynamic efhciency 
(ratio of the cooling-power to the absorbed power), the 
minimal temperature it allows and its applicability un¬ 
der diverse conditions. Here we examine these thermo¬ 
dynamic characteristics for a simple method we have 
previously introduced [IMl, whereby heat is trans¬ 
ferred from a reservoir by laser-assisted collisions to two- 
level atoms, and is subsequently emitted via atomic flu¬ 
orescence (Fig.[TJi). This method, hereafter referred to 
as laser-induced collisional redistribution (LICORE), is 
shown to be advantageous compared to sideband cooling 
BUI, regarding the ability to cool down a broadband 
heat reservoir. This ability is put to an experimental test 
for a hot bath of helium in collisional equilibrium with 
rubidium atoms and shown to be in good agreement with 
theory. 


The present work is essentially a thermodynamic anal¬ 
ysis of the cooling experiment performed in (lol - [l^ . Its 
aim is to compare theory and experiment in LICORE se¬ 
tups and focus on its main thermodynamic implications. 

Theory The LICORE scheme consists of an ensemble 
of two-level atoms, driven by the laser, and permanently 
coupled by elastic collisions to a hot bath, and by spon¬ 
taneous emission to a cold bath, the electromagnetic vac¬ 
uum, which is effectively at zero-temperature (Fig. [T^). 

For a hot bath at temperature T, the auto-correlation 
(coupling) spectrum, G//(w) = (t)^//(0))dt 

satisfies the detailed-balance (Kubo-Martin-Schwinger, 

KMS) condition [^: Gnioj) = Bnit) be¬ 

ing the bath operator in the interaction picture and fcg 
the Boltzmann constant. This detailed balance condition 
is akin to the Kennard-Stepanov ratio of absorption and 
emission rates (l6l - [l9l| . The spectral shape of Gnio}) is 
immaterial to the occurrence of cooling, but its peak and 
width affect the cooling rate, as shown below. 

The steady-state solution of the master equation (see 
Appendix A) defines an effective temperature Ttla, 


which is a measure of the stationary atomic level- 
populations ratio pee/Pgg, where Pee{gg) IS the excited- 
(ground-) state population of the two-level atom (TLA). 
The heat-flow direction between the atoms and the hot 
bath is determined by Ttla- For large and positive A, 
i.e. red detuning from the atomic resonance (A ^ 5), its 
Boltzmann factor satisfies 


Pee 

Pgg 


VpC '“Bl’ 
Fp + 7 



( 1 ) 


where the collision-induced pumping rate is Fp = 
Gl/(|A|). g is the laser-atom resonant coupling and 
7 is the spontaneous-emission rate. Hence, for A ^ g > 
0, the atoms are effectively colder (Ttla < T) than the 
hot bath. This is consistent with the expressions ob¬ 
tained from our general theory (l^ for the heat current 
flowing from the hot bath to a two-level atom (see SI, 
Eq. (S20) for the exact expression) 


Jh = h/^Tp - / _ aA X > 0. (2) 

7 -I- (^I + e j Fp 

Gonsequently, heat flows to the atoms, thereby cooling 
the hot bath, provided A > 0, for any T, 7 and Fp. 
The asymmetry between cooling and heating rates of the 
buffer gas in LICORE, as confirmed by the experimen¬ 
tal data (Fig. [2]), is expressed by 


-Jg(-A) _ ^ 

Jh{A) 


(3) 


This asymmetry is a consequence of two factors: 1) For 
|A| ^ g the coupling strength to the cold bath (elec¬ 
tromagnetic field) is independent of the sign of A. 2) 
The hot -bath (buffer gas) coupling to the TLA, GLf(A), 
satisfies the (KMS) detailed balance relation. The asym¬ 
metry in m conforms to the common intuition that it 
is easier to heat up than to cool down. LICORE can 
be understood as follows. Collisions broaden the atomic 
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Figure 1. (Color online) (a) Schematic representation of 
cooling by laser-induced collisional redistribution (LICORE). 
Left- the level scheme along with the corresponding laser de¬ 
tuning A, collisional pumping rate Fp, and the atomic spon¬ 
taneous decay rate 7 . Right- heat flow chart (TLA: two-level 
atom; EM: electromagnetic vacuum; H bath: hot bath; C 
bath: cold bath), (b) Sideband cooling compared to LICORE. 
Top: Efficiency as function of the band-width 'y^oU scaled by 
fl. While sideband cooling (solid red line) is highly efficient for 
resolved sidebands {Q S> 'Ycoii), LICORE (dashed blue line) 
is much more efficient for unresolved sidebands (fl <C 'Jcoii)- 
In sideband cooling fl is the trapping frequency, while in 
LICORE fl is the coupling Rabi frequency, taken to have 
the same value. The resolved sidebands spectrum conforms 
to the Mollow triplet HQ. 


levels, allowing the atom to absorb a photon below reso¬ 
nance frequency (A = wq — > 0). Energy flow from the 

hot bath compensates for the energy mismatch between 
absorption at v and spontaneous emission at uq to the 
electromagnetic vacuum (cold bath). The hot bath tem¬ 
perature is thereby reduced. In a perfectly isolated setup, 
this method would cool the hot bath down to an arbitrar¬ 
ily low temperature. As the hot bath temperature goes 
down, the cooling power in ([^ is reduced but remains 
always positive. This is a salient advantage compared 
to both Doppler cooling (^ . [2l| and sideband cooling 
[MQ that have fundamental minimal cooling tem¬ 
peratures (see Table 1). The LICORE cooling efhciency 
is a product of two factors, corresponding to a two-step 
process: Eirst, a photon is absorbed by the atom with 


efficiency Pl being the incident laser power. Next, 
the absorbed photon drives the cooling with the thermo¬ 
dynamic efficiency Thus the total efficiency in the 

weak laser limit (A ^ g, using ([2])) 


„ = = ^ _ (4) 

Pl Pabs Dl ^ Tp ’ 

the absorbed power being (see SI) 


hA 

^abs = Ep j RA \ ■ (^) 

7 -I- -I- e ^ j Tp 

The extracted energy from the hot bath follows the de¬ 
tuning, and correspondingly the cooling power increases, 
for moderate detuning in the linear regime (Eig. [2]). Yet, 
for large detuning laser absorption decreases, thereby 
reducing the cooling power. The highest efficiency is 
achieved if all photons are absorbed, (Pl = Pahs = 

), with unity efficiency being reached if the 
energy shift from the redistribution becomes as large 
as the entire cooling light photon energy, hv, yielding 
Jh = Pabs- The efficiency of such a heat distributor 
can in principle be impressively high (Eig. [T})). On the 
other hand, optical transitions in atoms do not allow for 
reasonable absorption rates at detunings of order of the 
absolute photon energy, and a more typical value for the 
high pressure buffer gas broadened ensemble is a detun¬ 
ing of order of the thermal energy, A Ri T, see Fig. 2, 
corresponding to A « 10 THz at 500 K gas temperature. 
This yields a typical efficiency Jn/Pahs ~ TI^o- 

To assess the merits of the LICORE method, we com¬ 
pare it (Table 1) to existing laser cooling methods: (i) 
Doppler cooling, which is based on the momentum change 
of an atom caused by the absorption and reemission of a 
photon (^; (ii) Sideband cooling which takes advantage 
of the sidebands created by the interplay between inter¬ 
nal and external degrees of freedom of a species Q to 
reduce its trapped-state occupation by tuning a laser to 
the lower sideband. The latter method requires (Fig.[Tl 3 ) 
the sidebands to be resolved, which amounts to strong 
binding to the trap D ^ 7 , otherwise heating processes 
will compete with the cooling. We can draw an analogy 
between the trapping frequency Q in sideband cooling 
and the coupling Rabi Frequency fl = y/Ag'^ -|- A^ in 
LICORE since in both processes the modes to be cooled 
have frequency O. Yet, we note that in LICORE the 
suppression of heating is achieved for any ratio of 7 /O 
provided there is weak eoupling to the laser, g A, 
whence O ~ A -|- , since the heating probability 

is then proportional to 

The thermodynamic efficiency bound determines the 
fraction of the absorbed energy used to cool down the 
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Table 1. Comparison of performance 



0 < 7 

D > 7 

Maximal thermodynamic efficiency ^ j 

Doppler Cooling f^. 



(m atom mass) <C 1 

Sideband cooling [3-[l, 

[3 

^B^min 1 1 

nn ~ T7rrT+774sTT ^ ^ 

7/40 j 

^B^min 1 1 

hO. ' , /'i+^2/i6sj2\ 

A/160^ ) 

Resolved bands — = —< 1 

Unresolved bands — = —R 7 T <C 1 

UJf\ — i 2 

Laser-induced colli- 

sional redistribution 

(LICORE) 

^R^min , 1 1 

''41og(f)^^ 

^B^min . 1 1 

Unresolved bands R A < 1 

V ~ 



Figure 2. (Color online) Experimental results (red dashed 
line) compared to theoretical prediction (blue solid line) of 
the total cooling power as a function of the laser detuning 
A. The y-axis indicates the cooling power. For red detun¬ 
ing, Jh corresponds to buffer gas cooling (purple, dark area), 
while blue detuning results in negative Jh , namely, buffer gas 
heating (orange, light area). 

bath. The Doppler cooling efficiency is low, because it 
scales with the recoil parameter. Sideband cooling and 
LICORE cooling may achieve similar efficiency for an ex¬ 
tremely narrowband oscillator bath, the scenario foreseen 
by sideband cooling, but LICORE may be highly efficient 
for a broadband bath, while sideband cooling efficiency 
will then be very low (Fig. [ 1 ) 3 ). 

Experiment The cooling setup of a high pressure cell 
with optical access contains a mixture of atomic rubidium 
vapor, with number density on the order of 10^®cm“^, 
and a noble buffer gas that serves as the hot bath, here 
helium gas with a number density of 1.5T0^^cm“^ to 
6T0^^cm“^. The cell temperature is kept around 500K. 

The confocal optical setup shown in Fig. [3] is designed 
to characterize the spectral shift of the emitted fluores¬ 
cence relative to the incident laser frequency. The col¬ 
lected atomic fluorescence radiation was analyzed with an 
optical spectrometer. Typical experimental data for an 
excitation laser frequency of 360 THz (401 THz), which is 
red (blue) detuned with respect to the D-liiies resonances 


Figure 3. Experimental setup for confocal spectroscopy on 
high pressure helium-rubidium gas mixtures. 


of the rubidium atoms, is shown in Fig. IT] The data 
displays typical pressure broadened fluorescence spectra, 
with the atomic rubidium Dl- and D2-lines clearly visi¬ 
ble, along with residual peaks of the exciting laser light. 
The figure clearly shows the redistribution of the fluores¬ 
cence frequency to the center of the D-lines for both red 
and blue detuning respectively. 

The cooling experiment is performed with a beam 
power of Pl = 2.4 W on a 1 mm beam diameter. On res¬ 
onance, the measured optical density in the 1 cm long gas 
cell, as derived from the transmitted laser power, reaches 
its maximum value of ca. 5. A more detailed description 
of the experimental setu p a nd the characterization of the 
system can be found in [ll|. 

In the present analysis we use experimental data for he¬ 
lium [l3| instead of argon [l^ as a buffer gas, as the ther¬ 
mal properties of helium, i.e. thermal conductivity and 
heat capacity differ both by an order of magnitude from 
the respective values in argon. Our theoretical model 
shows for both cases (helium and argon) excellent agree¬ 
ment with the experimental data, which demonstrates 
that the physical mechanism of cooling here is largely in¬ 
dependent of those quantities. Helium, due to its higher 
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thermal conductivity may be better suited than other no¬ 
ble gases to couple the cooling scheme to a cooling load. 

Comparison of theory and experiment We have 
measured the center of the fluorescence line and a{iy), 
the photon absorption probability. Upon identifying 
Pabs = Pha{v), Jh = PLa{v)^. We relate the theo¬ 
retical cooling power, Eqs @ - dSl), to its experimental 
counterpart. In our comparison, we account for laser ab¬ 
sorption by many rubidium atoms in the cell. 

The total cooling power can be calculated as: JH,tot = 
^ e~°‘^ JHiz)dz, where Na/L is the linear rubidium 
atomic density, and the laser power is assumed to expo¬ 
nentially decay with z along the cell due to absorption. 
For a typical cooling laser frequency of 365 THz, the mea¬ 
sured total absorption in the cell is 67%, corresponding 
to an absorption length /abs = l/a — 9 mm. 

A quantitative comparison with theory requires an 
ansatz for the coupling (autocorrelation) spectrum to 
the buffer gas bath, the simplest being a constant func¬ 
tion: Gh{uj > 0) = G'^r(O). This ansatz means that 
the spectrum of collisions that couple the buffer gas and 
the rubidium atoms is flat within the absorption line- 
shape of the laser, i.e. collisional effects are indepen¬ 
dent of A. The factor Gh{0) is determined from exper¬ 
imental data. We are aware that a more refined theory 
should account for the quasimolecular potential curves 
of rubidium-helium atom pairs [2^. Clearly, heat leak¬ 
age (non-isolation) or molecule formation in the collisions 
may stop the cooling. Experimentally, much lower tem¬ 
peratures are attainable provided that the cooling region 
is well isolated from heat leakage. 

Figure [2] compares the theoretical (blue solid line) and 
measured (red dashed line) cooling power. As predicted, 
the buffer gas (hot bath) is cooled down or heated up 
for red-detuned or blue-detuned laser, respectively. The 
difference between the measured and calculated curves 
reflects the difficulty to precisely determine the number 
of atoms to which the cooling power is proportional. Nev¬ 
ertheless, the main experimental features, i.e. the cooling 
and heating range and their asymmetry, are in very good 
agreement with theory. 

Minimal temperature predictions In principle 
much lower temperatures are attainable by LICORE 
than those observed by us, provided the technical issues 
discussed above are experimentally addressed. We again 
emphasize that the present theory does not account for 
quasimolecular character of the alkali-noble gas system. 
LICORE cooling may reach lower temperatures than 
sideband or Doppler cooling, including ultralow tempera¬ 
tures ( ^ 1 ) for both narrow and broadband spec¬ 

tra. Namely, D ~ A, can be either narrow or broad com¬ 
pared to 7 and still yield very low temperatures, while 
sideband cooling can reach such temperature only for re¬ 
solved sidebands (Table 1) O ^ 7 . Our theory predicts 
that for a given value of 0 at large detunings LICORE 



Frequency [THz] 

Figure 4. (Color online) Fluorescence spectra for rubidium 
vapor at 530 K with 200 bar helium buffer gas for excitation 
frequency 360THz (red) and 401 THz (blue). Arrows indicate 
the direction of the mean frequency shift of the scattered pho¬ 
tons by the redistribution process towards the D-lines center. 


may reach lower temperature than Doppler cooling (Ta¬ 
ble 1 ). 

In a perfectly isolated cell cooling stops only when 
Jh = 0 at the minimum buffer gas temperature obtain¬ 
able from the exact expression (SI) 


__n_ S+e -h 

e =- — 

S+ -h S-e~^ 


( 6 ) 


where the following short-hand notation is used: iy± = 
iz ± D, ^ For ^ weak laser and 

positive detuning (0 < 5 ^ A <C rz), iz+ Ri Wo and 
1 > 




the value 


3 > e , the minimal temperature attains 
-(see Table I). However, the 

41og(f) 

LICORE cooling efficiency and rate are very low at such 
large detunings, as seen from Eqs. (4) and (5) (Fig. lb). 

Discussion The main advantage of LICORE is that 
both the efficiency and the minimum temperature are 
adjustable by tunable laser parameters, A and g, while 
for Doppler or sideband cooling they depend on fixed 
parameters, spontaneous emission rate 7 and trapping 
frequency D. Furthermore, LICORE cooling does not 
impose any restriction on either the hot- or cold- bath 
spectra, and is linear in the laser power. It thus stands in 
stark contrast to refrigeration models 23|, [2^ where not 
only the modulation rate has to be above a critical rate, 
but also non-overlapping (resolved) hot and cold bath 
spectra are required (23|. However, the present scheme is 
a heat distributor and not a refrigerator, as it cools down 
the hot bath (rather than the cold bath) at the expense of 
the laser power. Without the laser irradiation, the effect 
of collisions with the bath would only be to decohere the 
atoms, rather than change the bath temperature. 
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Our combined experimental and theoretical results 
offer new fundamental and applied insights into the 
problem of laser-induced cooling of a heat bath. The fun¬ 
damental insight is that it is possible to run a laser-driven 
cooler that is not constrained by the thermodynamic 
bounds governing a refrigerator or a sideband cooler, 
and therefore allows cooling a broadband bath down to 
much lower temperature, and uniquely high efficiency. 
The applied insight is that a new technological avenue 
may be opened for cooling diverse gaseous or condensed 
media by laser-driven two-level atoms. Namely, the 
present method can in principle bypass the limitations 
of currently used Doppler or sideband cooling methods 
in terms of performance and versatility: the technique 
does not resort to auxiliary atomic levels or to resolvable 
sidebands', it allows for arbitrary (particularly, broad¬ 
band) bath spectra and temperatures and only requires 
a red detuning of the cooling laser. 
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APPENDIX 
A. Master Equation 
The total Hamiltonian is given by 

H = Hs{t)+ Y. \(Hsb)j + {Hb),1 (Al) 

3=H,C 

where the {Hb),j is the j-bath free Hamiltonian. Here 
the laser-driven system Hamiltonian is: 

Hs{t) = ^ojoaz + (A2) 

where ujq is the (resonance) frequency of the two-level 
atom, az and cr± — ax ± iay are the appropriate spin- 
1/2 operators, v is the laser frequency and g is the cou¬ 
pling strength between the laser and the two-level atom. 
The laser detuning is A = wq — The two-level atom 
coupling to the hot bath (H), via elastic collisions that 
do not change the two-level atom level populations, is 
described by 

{Hsb)h = haz ® Bh, (A3) 

where Bh is the buffer gas bath operator. The coupling 
to the cold (G) bath (electromagnetic vacuum) via spon¬ 
taneous emission is given by 

{Hsb)c = ® Be, (A4) 

Be being the electromagnetic bath operator. 

At steady-state, we may replace (|A2p of the main text 
by the averaged Hamiltonian 

Hs = (A5) 
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where A = loq — expresses the laser detuning. Then 
the coupling operator to the jth bath, Sj, is decomposed 
into a Fourier series 


Pabs = —Jc{t) — Jnit) 


(A14) 


S,it) = U\t)S,Uit) = 

^-i{Q+qv)t 

9=0.±l,±2...ez {w} 


^ {j=H,C) (A6) 


The master equation has the following form with the 
Lindblad generator decomposed in a Floquet (harmonic) 
series 

ps = Cps, (A7) 

i=ff,C(j=0.±l,±2... bJ 


^qaP = 

\{Gj{Cd + qv){[Sj^q{io)p,S]^^{Cd)] + ,(w),pS']_,j(<:;)]) + 

Gj(-w - (l'^){[S]^q{(P)p,Sj^q{u)] + [S']_,j(a;),pS'j,,(w)])|. 

(A8) 

Explicitly, in terms of fl = \/Ag'^ + A^, the laser- 
induced Rabi frequency, these coupling operators have 
the from 


Sc{v -n)^ 


A-VL / 0 1 
2VL loo 


Q ! \ 5/10 

= n 0 -1 


Sc{v + n) = 


A + Vt / 0 0 

2n I 1 0 


(A9) 


The population probabilities for the two-level atoms 
satisfy the following Markovian master equation (in the 
interaction picture) 

do ^ 1^1 

-^ =-(rp-l- 7 )pee + rpe“'^Pgg A >0 (A 15 ) 

and for negative detuning 

do ^1^1 

= -TpPee + {Tpe~'^ + 'y)pgg A < 0 (A 16 ) 


At steady-state Eq. (|A15ll yields Eq. (1) of the main 
text. 


B. Heating regime 

It is instructive to consider, for comparison sake, the 
heating regime obtained for A < 0. In this regime Ttla 
is determined by 


>i|A| 


e ^‘b'^tla = 


_ Pe 


BIAI 
Fi,e '‘bt 


■ T ^ 1^1 

— > e . 


Pgg 


(A17) 


In this case the atoms are hotter (Ttla > T) than the 
hot bath, thus the heat flows from the atoms, heating up 
the bath. In this case, the absorbed photons have more 
energy than the spontaneously emitted photons (wq > b') 
and the extra energy is absorbed by the bath, heating it 
up. The heat current is then 


By contrast 

J -1 


^ 1-10 


(AlO) 

The heat currents flowing out of the j-th bath, which 
are the main observables of interest, obey the same har¬ 
monic (Eloquet) expansion. They are given by 


jH = -n\A\rp 


/ ft|A| X 

7 -b (1 -b e ’‘B^ j 


< 0 . 


r„ 


(A18) 


The physical process is similar to that in the cooling case 
except that the Boltzmann factor in (2) of the main text 
is here absent, causing an asymmetry of the cooling and 
heating rates; 


Mt)= E (All) 

q— —1,0,1 CO 


Jqo(t) = -*B7(;Tr[(£G(i)p(t))inpA(t)]^ (A12) 


exT)/-lAh2iAH.l 

UjkBTi / 

PqCb ~ - 


Trexp{-(|gM}' 


(A13) 


C. Minimum temperature 

The most general expression for the heat current flow¬ 
ing from the hot bath is [l^ 

Jh = 


NMT^ 


e '^H(S+ + 6-e ^c)-(5+e ^^-bi5-) 

6- -b -b i5+ -b ^ 

(A19) 


The absorbed power is 


Under the conditions specified for Eq. (6) 
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— 

e 



(A20) 

D. Experimental setup 

The setup is shown in the main text (Fig. 3). It is a 
high presure cell with optical access. It contains a mix¬ 


ture of atomic rubidium vapor (whose density is of the 
order of 10^®cm“^) with helium gas (whose density goes 
from I.5T0^^cm“^ to 6T0^^cm“^). The cell walls are 
kept at around 500K. They are an extra heating source 
and its thermal contact with the buffer gas precludes 
the achievement of lower temperatures. Light from a 
titanium-sapphire laser is used to address the rubidium 
DI- and D2- resonances (near 377 THz and 384 THz, re¬ 
spectively) . Due to the high pressure of the helium buffer 
gas, the resonances are broadened to linewidths of a few 
nanometer. Nevertheless, we may treat the rubidium 
atoms as laser-driven two-level atoms at the chosen (D- 
line) transition. 



